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Abstract: Between September 2006 and February 2007, the galactic binary LS I +61 303 was monitored 
in the TeV band with the VERITAS array of imaging Cherenkov telescopes. These observations confirm 
LS I +61 303 as a variable TeV gamma-ray source, with emission peaking between orbital phase 0.6 and 
0.7. During this observational period, monitoring in the X-ray regime was also carried out using both the 
RXTE and Swift detectors, which offered complementary coverage of the source. Outbursts in the 0.2-10 
keV band were observed by both satellites at close to the same orbital phase as the TeV peak during the 
2 orbital cycles covered simultaneously in both bands. While this source has been extensively studied in 
the X-ray band in the past, this is the first observational campaign to utilize contemporaneous X-ray and 
TeV data on LS I +61 303. 



Introduction 

First associated with the COS-B source 2CG 
135+01 [1] in 1978, LS I +61 303 has been a source 
of considerable interest for the last 3 decades due 
to its peculiar behavior in radio and X-rays. Lo- 
cated at ~2 kpc distance [2], LS I +61 303 is a pair- 
ing of a Be star with a neutron star or a black hole 
which completes an orbital transit every 26.496 
days [3]. Flaring behavior has been historically de- 
tected in radio[4], X-rays[5][6], and most recently 
TeV gamma rays [7]. X-ray campaigns show flar- 
ing activity occuring regularly, every orbital cycle 
around orbital phase 0.4^0.7 [5] [6] [8] [9]. The 
most interesting feature of the X-ray flares is that 
they consistently appear to precede the radio out- 
bursts which occur several days later LS I +61 
303 has also been associated with the EGRET GeV 
source 3EG J0241[10] which showed evidence for 
outbursts both near perisastron passage (~0.23[3]) 
and later in the orbit near phase 0.5 [11]. 

Traditionally, there have been two classes of 
models surrounding variable high-energy emission 
from this source. The first [12] [13] assumes that 
the broadband emission is generated by accre- 



tion around the neutron star as it encounters vary- 
ing stellar wind densities in its orbit (microquasar 
models). Higher accretion rates cause a radio jet 
to form which can then upscatter stellar photons 
to TeV energies. Observations resolving compact 
one-sided radio emission from the system were 
taken as unambiguous evidence for the presence of 
a relativistic jet and thus proof of the microquasar 
model [13]. 

Alternatively, in the models of [14][15][16], the 
system is treated as a binary pulsar with particle 
acceleration taking place in the shock formed be- 
tween the pulsar and stellar winds. These pul- 
sar models interpret the resolved radio emission in 
[13] as coming from the cometary tail of the pulsar 
wind as it is blown about by the stellar wind. This 
assertion is lent further credence by the more in- 
depth radio imaging of this cometary tail detailed 
in [17]. Although the binary pulsar model for LS I 
+61 303 seems much more likely in light of these 
recent observations, there are still unresolved is- 
sues (see [18]) in both models which need further 
multiwavelength observations to investigate. 
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Figure 1: VERITAS, RXTE, and Swift monitoring of LS 1+61 303 over seven orbital periods. Peaks can 
be seen in aU three data sets near dates corresponding to (^=0.7 (also see figure 2). 



Observations 

VERITAS TeV Observations 

Sensitive in the energy range 100 GeV to 30 TeV, 
VERITAS is an newly commissioned array of four 
12-m imaging Cherenkov telescopes which has be- 
gun its full observational program as of January 
2007 (see G. Maier in these proceedings). Dur- 
ing the Fall of 2006, the initial 2-telescope stereo 
array was in operation, with a third telescope be- 
ing added for observations from January onward. 
TeV rates quoted henceforth are in terms of Crab 
Nebula fluxes taken from contemporaneous obser- 
vations (see O.CeUk in these proceedings). 

From September 2006 until February 2007, LS I 
+61 303 was observed for a total of 43.3 hours 
covering orbital phases 0.2— >0.9. These observa- 
tions covered 5 orbital cycles with reliable detec- 
tions resulting in every orbital period with the ex- 
ception of the February dark run (see G. Maier 
in these conference proceedings). However, the 
data taken in February were from a shorter ex- 
posure and are consistent with the behavior seen 
in the other months. During the "OFF" phases, 



i.e., 0.2^0.5 and 0.8^0.9, the source is quies- 
cent in the TeV range, with no reliable detection 
resulting. During flie "ON" phase, 0.5^0.8, VHE 
gamma rays are detected with emission ranging be- 
tween 4% and 13% of the Crab Nebula flux. VER- 
ITAS observations showed the TeV flux peak cen- 
tered at around 0.65^0.74, which is later in phase 
than the TeV maximum observed by MAGIC [7] at 
0.5^0.6. 

RXTE Observations 

Monitoring in the 2-10 keV band was carried out 
by RXTE [19] for nine 1-ks observations from 
2006/10/15 to 2006/10/31. Spectral fitting was car- 
ried out with XSPEC 5.12.3.1 assuming an ab- 
sorbed power law model. Fluxes shown in fig- 
ures 1 and 2 come from integrating this model 
over the 2-10 keV range and are scaled to emis- 
sion seen by RXTE from the Crab Nebula [20]. 
The light curve produced from this fitting shows a 
clear flare evolving from around orbital phase 0.45, 
peaking at around 0.7 and persisting until the end 
of the data set at 0.8. This flare was particularly 
strong, with emission peaking at 27xl0~^^ergs 
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Figure 2: Comparisons between contemporaneous TeV and X-ray data taken with VERITAS, RXTE, and 
Swift. TeV detections/upper limits are 3cr/99.9% Helene upper limits. All three observation sets show 
elevated emission between phases 0.5 and 0.8. 



s~^cm~^ which is larger than, but consistent with, 
past RXTE observations of this source[5]. The 
spectra fit to each data point show no strong evi- 
dence for spectral variation, however, the relatively 
short exposure time and reduced instrument oper- 
ation (i.e. only one PCU on for some nights) must 
be taken into account. A more detailed descrip- 
tion of these observations wiU be provided in an 
upcoming pubUcation. 

Swift Observations 

From September 2006 to October 2006, and again 
in November 2006-January 2007, the XRT instru- 
ment aboard the Swift satellite [21] was used to 
monitor LS 1 +61 303 in the 0.2-10 keV band. 
Reduction and analysis procedure can be found in 
Holder, Falcone and Morris (this conference). The 
XRT count rate to flux conversion factor used for 
this data was 5.7x10"^^ (erg/cm^/s)/(c/s). This 
conversion factor was obtained using the mean 
spectral fit parameters for these data. The system- 
atic error associated with spectral variation is not 
included in the error bars. This data set overlapped 



only once with VERITAS TeV data but fortunately 
it covered a region where VERITAS saw a clear de- 
tection of the source. Similar to the RXTE obser- 
vations, the Swift data set shows a clear elevation 
of emission occuring between orbital phases 0.5 
and 0.8. The data set also shows a secondary el- 
evated emission occurring between orbital phases 
0.0 and 0.1. 

Results 

TeV Variability 

It can be seen from figures 1 and 2 that the phase of 
the peak emission in the VERITAS data is always 
near the phase region 0.65^0.7. The VERITAS 
data have a maximum flux between 0.65 and 0.75, 
which, if compared to the location of the reported 
MAGIC TeV maximum of 0.5 to 0.6, appears to be 
shghtly offset in phase. However, a full statistical 
analysis of the light curves has yet to be carried out. 
Also, both the MAGIC and VERITAS light curves 
do not evenly sample the entire orbital phase. More 
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long-term, evenly sampled observations are needed 
to investigate the nature of the TeV variabiUty from 
the LSI +61 303. 

Correlation Between Data Sets 

Shown in figure 2 is a comparison between the 
TeV and X-ray data binned by orbital phase. It ap- 
pears that there may be a correlation between the 
rise time of TeV and X-ray emission in the sys- 
tem, although a full statistical analysis has yet to 
be carried out. Further analysis of the data as well 
as further long-term simultaneous observations are 
necessary to investigate this correlation further. 

Discussion 

With the prehminary analysis presented here we 
are unable to make any firm statements constrain- 
ing either the microquasar or the binary pulsar 
models. It is noted that the X-ray and TeV emis- 
sion appear to be correlated in time. Further anal- 
ysis of this data set will be performed to study 
the nature of a possible correlation between the 
collected data from RXTE, Swift, and VERITAS. 
This analysis will be presented in an upcoming 
pubhcation. 
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